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Abstract 


We  have  utilized  the  ligand  binding  properties  of  a  buried  cavity  created  in  the  interior  of  a 
protein  to  obtain  direct  information  about  the  variation  in  the  strength  of  CH-0  interactions 
between  the  ligand  and  protein.  Our  study  shows  that  the  strength  of  CH-0  interactions  can  be 
modulated  by  over  1  kcal/mol  by  changes  in  the  C-H  band  polarity.  Consequently,  several  such 
interactions  may  play  a  significant  role  in  the  stability  of  macromolecular  structures. 
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1.  Introduction 


While  hydrogen  bonds  of  the  type  XH-Y  (X,  Y  =  N  and/or  O)  are  essentially  found  in  all 
macromolecular  structures,  the  role  of  CH-0  hydrogen  bonds  in  such  systems  is  uncertain.  Recent 
reports  [1-5]  have  suggested  that  a  significant  number  of  the  CH-0  contacts  observed  in  proteins, 
RNA,  and  carbohydrates  represent  cohesive  interactions,  which  may,  in  sufficient  numbers,  be 
significant  to  structure,  stability,  and  function.  The  existence  of  short  intermolecular  CH-0 
interactions  is  well-established  in  many  small-molecule  crystals  [6-9].  In  cases  where  the  C-H-0 
angle  is  approximately  linear  and  the  C-0  bond  distance  is  less  than  the  combined  van  der  Waals 
radii  (3.3  A),  these  interactions  have  been  labeled  as  hydrogen  bonding.  Theoretical  calculations 
have  estimated  the  CH-0  bonding  interactions  to  be  worth  1-2  kcal/mol  [2, 10,  11].  However, 
many  computational  approaches  utilizing  force-field  parameterizations  do  not  explicitly  consider 
such  interactions  as  attractive,  and  CH-O  hydrogen  bonds  are  not  generally  considered  in  the 
analysis  of  structures.  The  strength  of  such  interactions  should  vary  considerably  depending  on  the 
C— H  bond  polarity  and  may  thus  be  different  for  the  C0H-0  interactions  observed  in  some  protein 
P-sheets  than  those  for  others  involving  the  more  acidic  C€H  protons  of  histidine  [1].  While  there 
is  a  correlation  of  the  C-O  distance  with  the  C— H  bond  acidity  for  a  number  of  organic  compounds 
[6, 12-14],  nothing  is  known  experimentally  about  the  strength  of  these  interactions  and  how  this 
strength  varies  with  the  polarity  of  the  proton  donor. 


2.  Background 


We  have  utilized  the  ligand  binding  properties  of  a  buried  cavity  created  in  the  interior  of  a 
protein  to  obtain  direct  information  about  the  variation  in  the  strength  of  CH-O  interactions  between 
ligand  and  protein.  The  cavity  created  by  the  W 191 G  mutation  of  cytochrome  c  peroxidase  (CCP) 
has  been  shown  to  bind  a  number  of  cationic  heterocyclic  compounds  [15-17].  Two  such 
compounds,  2,3,4-trimethylthiazole  (234TMT)  and  3,4,5-trimethylthiazole  (345TMT)  bind  to  the 
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W191G  cavity  in  essentially  isosteric  conformations  (Figure  1).*  Crystal  structures  of  the  protein 
soaked  in  solutions  containing  these  compounds  show  clear  evidence  for  binding  in  the  omit  electron 
density  maps.  No  significant  differences  were  observed  in  the  structure  of  the  protein  between  the 
ligand-free  and  either  of  the  ligand-bound  states.  Placement  of  a  ligand  model  within  the  omit 
density  implies  a  close  contact  between  the  C5  ring  proton  of  234TMT  and  one  carboxylate  oxygen 
of  Asp-235.  An  analogous  interaction  involving  the  C2  ring  proton  is  observed  for  345TMT. 
Asp-235  is  observed  to  hydrogen-bond  to  other  cavity-binding  ligands  and  helps  determine  the  cation 
bonding  specificity  ofthis  cavity  [15, 17, 18].  Thus,  the  absence  of  standard  hydrogen-bond  donors 
in  234TMT  and  345TMT  evidently  results  in  the  substitution  of  weaker  alternative  interactions  that 
fulfill  a  similar  role. 

The  geometries  of  the  interactions  implied  by  the  crystal  structures  are  consistent  with  a  CH-0 
hydrogen  bond  [8].  The  C-O  distances  between  the  ligands  and  the  Asp-235  carboxylate  oxygen 
were  2.95  and  2.85  A  for  234TMT  and  345TMT,  respectively,  shorter  than  the  combined  van  der 
Waals  radii  of  C  and  O  (3.3  A).  The  C-H-0  angles  estimated  from  placing  the  geometrically 
optimized  ligandst  into  the  electron  density  were  approximately  150°  and  140°  for  234TMT  and 
345TMT,  respectively  (Figure  2).  The  deviations  of  these  angles  from  an  ideal  linear  CH-  0 
hydrogen  bond  are  not  outside  the  range  (130-170°)  of  those  observed  in  organic  crystals  [8]. 
Although  the  two  compounds  differ  structurally  only  by  the  interchange  of  the  N3  and  C4  atoms  in 
the  thiazole  ring,  this  difference  serves  to  increase  the  acidity  of  the  C2  proton  in  345TMT 


*  Crystal  structures  of  W191G  soaked  in  50-mM  234TMT  or  345TMT  were  determined  as  previously  described  [17]. 
For  each  ligand,  the  orientation  was  unambiguous  due  to  density  observed  for  the  methyl  substituents  and  for  the  sulfur 
atom  at  the  higher  contour  level.  The  maps  were  constructed  by  merging  for  the  soaked  crystal  with  from 
models  of  the  W191G  empty  cavity.  The  structures  of  234TMT  and  345TMT  bound  to  W191G  have  been  submitted 
to  the  Brookhaven  database  (entries  lac4  and  lac8,  respectively).  Diffraction  statistics  for  234TMT  and  345TMT, 
respectively,  included  unit  cell  parameters  (103.76, 73.31, 44.64  A  and  106.03, 75.89, 51.13  A),  resolution  (2.0  and 
2.1  A),  Ilol(tv)  (18.2  and  12.4),  Ilcr,^  *6®  (2.00  and  1.76),  number  of  reflections  (18,219  and  19,056),  percent 
completeness  (80  and  87%),  and  RiyJD  (0.048  and  0.062). 

+  Ligand  partial  charges  were  calculated  as  Kollman  electrostatic  potential  (ESP)  charges  from  geometry-optimized 
structures  using  density  functional  theory.  The  calculations  were  performed  in  vacuo  as  previously  described  [19] 
using  the  program  Gaussian  94,  Becke3LYP  functional,  and  6-3 1G*  basis  set  ESP  partial  charges  were  used  to 
generate  an  ESP  grid  with  the  program  Delphi  (MSI)  using  a  uniform  dielectric  constant  of  1,  ionic  strength  of  0.0, 
and  full  Coulombic  boundary  conditions.  This  potential  grid  was  mapped  onto  the  solvent-accessible  surface  of  the 
molecule  calculated  with  a  1.4-A  probe  radius  as  a  color  spectrum  from  +135  kT/e  (red)  to  +200  kT/e  (blue). 
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Figure  1.  Electron  Density  for  (A)  234TMT  and  (B)  345TMT  Binding  to  the  W191G  Cavity. 

Shown  Are  Stereoviews  of  F0-Fc  Fourier  Omit  Maps  Contoured  at  +4  a  (White)  and 
+9  a  (Red)  Superimposed  on  a  Model  of  the  Ligand  That  Was  Placed  Into  the  Omit 
Density.  (See  Footnote  [*]  on  p.  2.) 

(p£a  *  17-19)  [20,  21]  relative  to  the  C5  proton  in  234TMT  (pATa  ~  25-30).  This  difference  in 
acidities  is  manifested  in  the  electrostatic  potential  calculated  for  the  two  compounds  and  mapped 
onto  the  solvent-accessible  surface  (Figure  2).  Since  the  CH-0  hydrogen  bond  is  primarily 
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Figure  2.  Electrostatic  Potential  Surfaces  for  (A)  234TMT  and  (B)  345TMT.  The  Color  Map 
Illustrates  the  Localized  Differences  Resulting  From  the  Increased  Polarity  of  the 
C2-H  Proton  for  345TMT  Relative  to  the  C5-H  Bond  of  234TMT.  (See  Footnote 
[t]  on  p.  2.)  At  the  Bottom  Are  the  Geometries  of  the  Interactions  With  Asp-235 
That  Are  Inferred  From  the  Placement  of  the  Geometry-Optimized  Ligand  Into  the 
Omit  Electron  Density. 


electrostatic  in  nature,  the  calculated  potentials  indicate  that  345TMT  should  form  a  significantly 
stronger  interaction  with  Asp-235. 


Thermodynamic  parameters  of  234TMT  and  345TMT  binding  to  W191G  were  determined  by 
isothermal  titration  calorimetry  (Table  1).  The  results  show  a  five-fold  decrease  in  Kd  for  345TMT 
relative  to  234TMT,  corresponding  to  1 .2  kcal/mol  of  additional  binding  energy  for  345TMT.  Due 


4 


Table  1.  Thermodynamic  Parameters  for  Trimethylthiazole  Binding  to  the  Buried  Cavity  of 
W191G* 


A H 

AS 

AG 

(mM) 

(kcal  mol'1) 

(cal  mol'1  K'1) 

(kcal  mol'1) 

234TMT 

1.50(0.16) 

-14(0.6) 

-33(2) 

-3.9  (0.1) 

345TMT 

0.20  (0.03) 

-15(0.7) 

-34(3) 

-5.1  (0.2) 

to  the  absence  of  structural  changes  in  the  protein,  the  similar  contacts  made  with  the  protein,  and 
the  fact  that  both  compounds  desolvate  the  cavity  to  the  same  degree,  the  observed  difference  in  the 
free  energy  of  binding  for  the  two  compounds  should  result  from  differences  in  the  electrostatic 
interactions  between  ligand  and  protein  and  from  differences  in  their  desolvation  energies.  However, 
estimates  of  ligand  desolvation  energies  were  very  similar  (Table  2),  indicating  that  the  observed 
difference  in  binding  free  energy  for  the  two  compounds  can  be  attributed  to  protein-ligand 
interactions. 


Table  2.  Calculated  Electrostatic  Contribution  to  the  Relative  Binding  Energy  of  234TMT 
and  345TMT  to  the  W191G  Cavity* 


-""a 

i 

> 

-avb 

-AV„ 

^^solv(  water) 

234TMT 

-56.88 

3.82 

-6.44 

-7.26 

-66.76 

-26.112 

345TMT 

-57.83 

2.98 

-5.26 

-7.13 

-67.24 

-26.096 

AAG345tmt_234Tmt 

-0.95 

-0.84 

1.18 

0.13 

-0.48 

-0.016 

NOTE:  Values  are  given  in  kcal/mol. 

The  electrostatic  contribution  to  the  ligand-protein  interaction  was  estimated  for  the  two  ligands 
bound  to  the  protein  cavity.  Partial  charges  were  calculated  for  each  ligand  (Figure  2)  to  account  for 


Dissociation  constants  and  enthalpies  were  measured  at  25°  C  in  100  mM  Bis-Tris  propane  pH  4.5  by  isothermal 
titration  calorimetry  (Microcal  MC2 ITC  calorimeter).  W191G  (0.2-0.4  mM)  was  titrated  with  5-pL  injections  of 
ligand  (2-8  mM)  equilibrated  in  the  same  buffer.  Error  estimates  for  Kd  and  A H  are  given  as  the  standard  deviation 
of  multiple  determinations,  and  those  for  AS  and  AG  were  obtained  by  propagation. 
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the  different  C-Hbond  polarities,  and  the  protein  dipoles/Langevin  dipoles  (PDLD)  method  [  1 9, 22]* 
was  used  to  account  for  effects  of  microscopic  atomic  polarizabilities  of  atoms  surrounding  the 
ligands.  The  more  favorable  electrostatic  interaction  of  345TMT  with  the  protein  relative  to  that  of 
234TMT  (Table  2)  arises  from  both  dipole  (AV^)  and  induced  dipole  (AVCo)  terms.  These 
interactions  are  partially  canceled  by  the  solvent  screening  effects  represented  in  the  Langevin  grid 
(A  VQL)  and  Bom  (AVB)  terms.  The  net  electrostatic  contribution  to  the  binding  free  energy  (AVtet) 
is  approximately  0.5  kcal/mol  of  additional  stabilization  of 345TMT  relative  to  234TMT .  This  value 
is  in  good  agreement  with  recent  ab  initio  quantum-mechanical  analysis  of  CH-0  bond  energies  [2] . 
Thus,  electrostatics  alone  form  a  significant  contribution  to  the  observed  difference  in  the  binding 
free  energy  of  1.2  kcal/mol. 

3.  Conclusion 

In  conclusion,  the  observation  of  untethered  ligand  binding  to  this  engineered  cavity  shows  that 
the  CH-O  interaction  is  made  by  choice  and  is  thus  a  stabilizing  interaction.  While  it  is  weak, 
summation  of  several  such  interactions  may  play  a  significant  role  in  the  stability  of  macromolecular 
structures.  The  strength  of  this  interaction  can  be  modulated  by  over  1  kcal/mol  by  changes  in  the 
C-H  bond  polarity.  Thus,  selective  inclusion  of  polar  CH-0  interactions  in  refinement  algorithms 
for  nuclear  magnetic  resonance  (NMR)  and  crystal  structures  may  be  justified. 


*  Electrostatic  calculations  were  performed  using  the  PDLD  method  with  the  program  POLARIS  [22].  AV10t  is  the  sum 
of  four  terms,  A V^,,  AVg,,,  AVgL,  and  AVB.  AV^  is  the  classical  electrostatic  interaction  of  the  ligand  charges  with 
all  of  the  charges  on  all  protein  atoms.  AV^  is  the  interaction  energy  from  the  dipole  moments  induced  by  the  electric 
field  as  a  result  of  atomic  polarizabilities.  AVCL  is  the  energy  of  interaction  of  the  ligand  with  the  field  defined  by  the 
solvent  dipoles.  A  VB  is  the  Bom  energy  of  interaction  with  bulk  solvent  beyond  the  radius  used  to  define  the  Langevin 
dipole  grid.  Values  of  A  V  have  the  sign  of  electrostatic  potential,  and,  thus,  a  more  positive  value  for  A  V  corresponds 
to  a  more  negative  AG.  POLARIS  calculations  were  performed  with  Asp-235  and  Arg-48  charged,  and  all  other  amino 
acid  residues  netural.  An  estimate  of  the  electrostatic  interaction  of  the  ligand  charges  with  the  solvent  reaction  field 
(A<7h*(MH)  for  the  two  compounds  was  obtained  with  a  continuum  dielectric  model  (Delphi,  MSI)  using  e,^  =  80 
and  =  2,  with  a  grid  spacing  of  0.25  A. 
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